Spontaneous persistent currents in magnetically ordered graphene ribbons 
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We present a new mechanism for dissipationless persistent charge current. Two dimensional 
topological insulators hold dissipationless spin currents in their edges so that, for a given spin 
orientation, a net charge current flows which is exactly compensated by the counter-flow of the 
opposite spin. Here we show that ferromagnetic order in the edge upgrades the spin currents into 
persistent charge currents, without applied fields. For that matter, we study an interacting graphene 
zigzag ribbon with spin-orbit coupling. We find three electronic phases with magnetic edges that 
carry currents reaching 0.4nA, comparable to persistent currents in metallic rings, for the small spin 
orbit coupling in graphene. One of the phases is a valley half-metal. 

PACS numbers: 



Ordered electronic phases can emerge in condensed 
matter with properties fundamentally different from 
those of the constituent atoms. Two main different sce- 
narios are known that result in the emergence of non- 
trivial electronic order. On one side, spontaneous sym- 
metry breaking driven by many-body interactions which 
accounts for the existence of the crystalline order in 
solids and the variety of ordered electronic phases they 
can present, like superconductivity and ferromagnetism 
[l[. On the other side, topological order, which accounts 
for the robust quantized properties of the electron gas 
in the Quantum Hall regimes, and, more recently, on 
the properties of the so called topological insulators 0- 
5]. Whereas Integer Quantum Hall state is driven by an 
external magnetic field, topological insulators are driven 
by spin orbit interaction. They are different from con- 
ventional insulators because of their conducting surface 
(or edge) spin states which can be either chiral or spin- 
filtered and are robustwith. respect to time-reversal sym- 
metric perturbations 



The experimental finding of topological insulators @,H| 
and their exotic surface (edge) states, motivates the gen- 
eral questions of whether and how electronic interactions 
could produce electronic phase transitions in the surface 
or edge states of a topological insulator, and what would 
be the consequences of such symmetry breaking. Here 
we address these questions in the case of a graphene 
ribbon with zigzag edges, a system that atracts enor- 
mous interest [3, in the context of spintronicsd, 

fl El and 
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0], magnetoelectronics[l 

valley tronics [l6| . When described with a one orbital 
tight-binding model, graphene zigzag ribbons are con- 
ducting because of two degenerate almost fiat bands as- 
sociated to states localized at the edges [9| . Importantly, 
the zigzag states preserve the valley character [111 . [l6j of 
two dimensional graphene. When considered separately, 
the effect of electronic repulsion and spin orbit interac- 
tions over the system is dramatic and has been widely 
studied. On one side, it was soon recognized that elec- 



tron Coulomb repulsion, in the Hubbard model, results 
in ferromagnetic order at the edges of the ribbon [lOj . In 
the ground state the edges are counter-polarized and ac- 
cording to both density functional calculations 12j and 
to mean field Hubbard model (2fjj. a gap opens so that 
the ribbon is an antiferromagnetic insulator. 

When Coulomb repulsion is neglected, spin orbit 
(SO) interaction opens a gap in the spectrum of bulk 
graphene 23J, |2J] which is accompained by the emergence 
of spin filtered edge states at the Fermi energy [j, 0]. 
These two features, a spin-orbit driven gap and the 
emergence of topologically robust edge states, are the 
main properties of Topological Insulators Thus, the 
Coulomb driven and the SO driven phase have very dif- 
ferent magnetic and conducting properties. Here we ad- 
dress the electronic properties of the ribbon when both 
Coulomb repulsion and spin orbit coupling are considered 
within a mean field Hubbard model[10|, |20( with Kane- 
Mele spin orbit coupling 0, We report two main find- 
ings. First, in the presence of SO coupling the states with 
counter-polarized ferromagnetic edges break valley sym- 
metry and, above a critical SO strength, the gap closes in 
one valley only, resulting in a valley half-metal. Second, 
in the presence of SO coupling ferromagnetic edges give 
rise to charge currents in the edges without an applied 
magnetic field. 

We consider the one orbital Hubbard model in a hon- 
eycomb zigzag ribbon, at half filling, with the addition 
of the Kane-Mele SO coupling[3J: 



(1) 



The hopping matrix Uj(a) in the first term accounts 
both for the standard first-neighbor spin-independent 
hopping t and the Kane-Mele-Haldane (KMH) second 
neighbor spin-orbit coupling [3, 0, HEf. The amplitude of 

the latter is given by the expression Ukm&z ■ (^di x , 
where d\ t 2 are unit vectors along the direction of the bond 
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FIG. 1: (a) Honeycomb zigzag ribon with N y = 6 atom 
rows. The second neighbour hopping vectors di are shown, 
(b) Low energy bands for the Spin Hall phase in N y = 12 
ribbon (U = 0,tKM = 0.01 eV). (c) Low energy bands for 
the same ribbon with counter-polarized ferromagnetic edges 
((7 = 3 eV and tKM = 0). (d) Contour map of the valence 
band wave function ^(^(j/)! 2 for spin up (left panel) and 
spin down (right panel), (e) Same than (d) for the magnetic 
ribbon. 



that connect atom i and j with their common first neigh- 
bor (see figure la), z is the unit vector normal to the rib- 
bon plane, and a = ±1 indexes the spin projection along 
z. For a flat ribbon, the SO term conmutes with a z . 
The second term in eq. ([1]) describes Coulomb repulsion 
between electrons in the Hubbard approximation. We 
treat it in a mean field approximation so that we end up 
with an effective single particle Hamiltonian where the 
electrons interact with a spin dependent potential that is 
calculated self-consistently, U J2i ( n i,T + n i,i-( n i,t))i 
so that a z is a good quantum number. Zigzag ribbons are 
defined by N y , the number of zigzag chains which yield 
a total of 27Y a atoms per unit cell in a one dimensional 
crystal (fig. la). Importantly, the top and bottom edges 
belong to the two different triangular sub-lattices that 
define the honeycomb lattice. For a given wavenumber 
k and spin a the mean field Hamiltonian has 2N y states 
^kau{y) with energy e au (k). 

We consider first the model in the two limit cases, 



U = 0, tuM > and txM = 0,U > for a ribbon 
with Ny = 12. The highest occupied and lowest empty 
energy bands of the U = 0, txM > case are shown in 
figure lb. The special spin filtered edge states are the 
linear bands crossing the Fermi energy (Ep) The wave 
function squared, \^k,a,v(y)\ 2 , of the valence band states 
are represented in figure Id for a =f and a =\. respec- 
tively. It is apparent that a =t {a =1) electrons can 
only be in the top (bottom) edge for positive (negative) 
velocity states. Spin f and I electrons can also be in 
the bottom edge, but with velocity opposite to that of 
the top edge. Thus, the ribbon is conducting, the edges 
are not spin-polarized, but the states at the Fermi en- 
ergy carry a net spin current. The electronic structure 
of the same ribbon, but now taking t^M = 0, U > 0, 
is radically different. Instead of the linear spin-filtered 
bands there is a gap, Ao so that the system is insulat- 
ing. The wave-functions of the valence band are shown 
in figure le. They do not show correlation between spin 
and velocity, but they are edge-sensitive: top (bottom) 
edge is ferromagnetic with spin up (down) majority. The 
solution with reversed spins is equally valid. 

In the light of figure 1 , the question of how these com- 
peting electronic phases merge when both the Coulomb 
repulsion and the SO coupling are present calls for an 
answer. In figure 2 we show the energy bands for three 
cases, all with U = t = 3eV and N y = 12. In the upper- 
panels we show two cases with magnetic edges with total 
spin zero, one with with txM = 0.014 and the second 
with txM = 0.03t. As shown in the inset, the magnetic 
moments are localized in the edges. In the small txM 
case, the system in an insulator, but the inter-edge gaps 
Ao are now valley dependent. As the SO coupling in- 
creases the gap in one of the valleys closes completely, 
yet the edges are magnetic. This phase is radically dif- 
ferent from the SO free case: a valley half-metal with 
magnetic edges antiferromagnetically oriented. In figure 
2d we show how the magnetic moment in the edge atoms, 
m = { n t)-( n 4.) is depleted as the strength of the spin orbit 
coupling is increased, reflecting the competition between 
the two terms in the Hamiltonian. Above a certain value 
of txni, the magnetic moment vanishes altogether and 
a Spin Hall insulator phase with conducting edge states 
identical to that with U — is obtained. The calculated 
phase diagram in the (tx m , U) plane showing the three 
different phases with zero total spin is shown in figure 
3a, for the case of Ny = 12. In figure 2c we show the 
conduction and valence band of the ferromagnetic phase 
with copolarized edges (see inset of fig. 2c). These bands 
preserve valley symmetry and intersect the Fermi energy, 
so that the edge states are conducting. 

The evolution of the valley symmetry breaking, re- 
flected by the different size of Ao(K) and Aq(K'), is 
shown in figure 3b. The valley symmetry breaking 
can be understood using perturbative arguments. Close 
to the Dirac point the Kane-Mele Hamiltonian can be 
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FIG. 2: (Color Online). Valence and conduction band of a 
N y — 12 ribbon with U — t — 3eV for three cases: (a) AF 
insulator with tuM = O.Olf, (b), AF valley half metal tfCM = 
0.034 and (c) Co-polarized ferromagnetic edges tKM = 0.03i. 
Insetsicalculated magnetic density along the ribbon cell, (d) 
Depletion of the edge magnetic moment as the SO coupling 
tKM is increased, for two different values oiU/t. 



approximated [3( by Vkm = ^V^tKM^zTz^zt where a z , 
t z and X z are the spin, valley and sub-lattice index re- 
spectively. Let us consider first the {txM — 0,U > 0) 
antiferromagnetic phase (figure lc,e) as starting point. 
The valence band is made of states with a z =f and sub- 
lattice X z = + (bottom edge) and states with a z =1 and 
sublattice X z = — (upper edge), so that in both cases 
the product a z \ z has the same sign. Now it is apparent 
that the expectation value of <j z t z \ z has opposite signs 
in opposite valleys, so that in one case the gap opens and 
in the other closes, as shown in figure 2a. 



b 1.5 



0.3- 





0.05 



0.1 



FIG. 3: (Color online), (a) Phase diagram for the AF in- 
sulating, AF valley half-metal (HM) and non-magnetic (NM) 
phases for the N y — 12 ribbon, (b) Evolution of the gaps 
Ao(K) and A(K') as a function of the SO coupling. 



In the non-interacting Kane-Mele model a non- 
equilibrium current / induces spin accumulation m in 
the edges 0. This can be quantified as follows. We as- 
sume that a population of non-equilibrium extra elec- 
trons occupies the positive velocity states, so that the 
Fermi energy is increased by 5Ep = eV. Using Landauer 
formula we have I — 2j^5Ep, half of which goes on the 
top edge. From the non-interacting conduction band dis- 
persion 6^ = hvpk, we obtain Skp and the corresponding 
change in density due to 5Ep, Sn = ^ — -^Skp. Since 
SEp is small, this extra density goes as spin f to one edge 
and spin .J, to the other. Thus, half of the extra charge 
goes to each edge fully spin-polarized, so that the edge 
magnetic moment reads 8m = jSq = j^Skpa. We can 
write the current in a given edge as 



evp 

Air m 

o 



(2) 



We now show that this picture survives in the inter- 
acting case in equilibrium. We find that in the ground 
state of the the magnetically ordered topological insula- 
tor phases shown in figure 2, charge currents flow on the 
edges, the total current accross a unit cell being null, in 
agreement with general theorems 26] . The current oper- 
ator is given by the sum of link currents associated with 
all the sites b connected to a by single particle hopping: 




(3) 



For a given eigenstate ^k a i> of the mean field Hamiltonian 
the current across the link ab reads 



lab [*] 
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(4) 



where <p ao — for a and b in the same cell, 4> a b = ika if 
a and b are in adjacent cells, and N is the number of cells 
in the crystal. The expectation value of this operator in 
the many-body ground state is obtained summing over all 
occupied bands: (I a b) = J2k.v,a f ( e ^( fc )) ^o6- Tne aver- 
age current is defined in the links of any pair of atoms 
connected by hopping in the one-body hamiltonian. In 
figure 4 we plot the ground state current map for the 
N y = 12 ribbon for the AF insulating phase (fig 2a) and 
the ferromagnetic conducting phase (fig. 2c). The AF 
valley-half metal (not shown) is very similar to the AF 
insulator. The three electronic phases described in figure 
2 present edge current of similar magnitude. Whereas 
in the FM case, currents flows in opposite directions in 
opposite edges, in the AF phases current runs parallel 
in the two edges. These results can be rationalized as if 
the magnetization plays the role of an external magnetic 
field. Thus, in the ferromagnetic case current flows is the 
same than in a Quantum Hall bar. In the antiferromag- 
netic phases, though, current flows parallel in the two 
edges, as if a magnetic field was pointing along opposite 
directions in the two edges. 
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FIG. 4: Current and magnetization maps for (a) the AF 
insulating N y = 12 ribbon and (b) the FM conducting one, 
both with tKM = 0.03i and U = l.Ot. The current between 
two atoms r a and fb is plotted as a vector along the line r& — r a , 
starting in the midpoint, (c), Edge charge current, in units of 
f tKM, as a function of the magnetic moment of the edge, for 
3 values of txM = 10 -2 , 10 _1 and lmeV, for the AF phases. 



pKM 



The magnitude of the top-edge current, normalized by 
as a function of the top edge magnetization m 
collapses for several values of tKM (figure 4c). For small 
m the curve is linear, I e d ge — —^itKMfn, in qualitative 
agreement with the analytical result of eq. ©, since 
vp = jtKMa/fo with 7 ~ 6.5. At larger m the edge 
current saturates to |/ e dge| — 0.4|tifM- If we take txn — 
lOfieV, close to the small values obtained by ab-initio 
calculations 27J, we obtain an edge current of ~ OAnA, 
well within reach of state of the art persistent current 
detection [28j|. 

In conclusion, we propose a new mechanism for persis- 
tent charge currents. It involves the edge states of a two 
dimensional topological insulator with spontaneous fer- 
romagnetic order induced by Coulomb interactions. We 
propose that this scenario occurs naturally in graphene 
zigzag ribbons. We find three new electronic phases in 
that system that combine ferromagnetic order and spon- 
taneous charge current flow, both localized in the zigzag 
edges. They arise from the interplay of Coulomb repul- 
sion and spin orbit coupling. When the ferromagnetic 
edges are counter-polarized, the valley symmetry is bro- 
ken and, above a critical strength of the spin orbit cou- 
pling, the system goes from an insulating to a valley- 
half-mctal phase. In the three phases, current flows as 



if there was a real magnetic field perpendicular to the 
sample along the direction of the magnetization edge. 
The question of weather our findings can be generalized 
to magnetically doped topological insulators will be ad- 
dressed in future work. 
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